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We show that by using a feedback loop it is possible to reduce the fluctuations in one quadrature 
of the vibrational degree of freedom of a trapped ion below the quantum limit. The stationary state 
is not a proper squeezed state, but rather a squashed state, since the uncertainty in the orthogonal 
quadrature, which is larger than the standard quantum limit, is unaffected by the feedback action. 
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j /^ I In recent years there has been an increasing interest on trapping phenomena and related cooling techniques [jj] . A 
number of recent theoretical and experimental papers have investigated the ability to coherently control or "engineer" 
atomic quantum states. Experiments on trapped ions, where the zero point of motion was closely approached trough 
laser cooling already showed the effects of nonclassical motion in the absorption spectrum ^||. More recent 
experiments report the generation of Fock, squeezed and Schrodinger cat states H,|j. These states appear to be of 
fundamental physical interest and possibly of use for sensitive detection of small forces Q. Moreover, the possibility 
to synthetize nonclassical motional states gave rise to new models in quantum computation j?J . 

In this paper we present a way to reach a stationary nonclassical motional state for a trapped particle which 
is able to give a significant uncertainty contraction in one phase-space direction. The scheme can be then applied to 
control the vibrational motion against the heating processes responsible for decoherence || . This could be important 
to obtain high fidelity in quantum logic operations || . 

The basic idea of the scheme is to realize an effective and continuous measurement of a vibrational quadrature for 
i-G | the trapped particle and then apply a feedback loop able to control, i.e. to reduce, its fluctuations even below the 
quantum limit. 

The fact that a feedback loop may reduce the fluctuations in one quadrature of an in- loop field without increasing the 
fluctuations in the other has been known for a long time, and has recently been called "squashing" |l(| as opposed to 
"squeezing" of a free field, in which the conjugate fluctuations are increased. In our scheme, the obtained stationary 
O^ 1 state results as a "squashed" state; however, the quadrature measurement increases the noise in the orthogonal 
quadrature well above the standard quantum limit and the squashing comes with respect to the state one has in 
presence of the measurement process and without the feedback action. In this way the uncertainty principle is not 
rS ' violated. 

The paper is organized as follows. In section II we show how to realize the indirect continuous measurement of 
a vibrational quadrature by coupling the trapped particle with a standing wave. In section III we shall introduce 
the feedback loop, in section IV we shall study the properties of the stationary state in the presence of feedback and 
section V is for concluding remarks. 



II. CONTINUOUS MONITORING OF ATOMIC MOTION 



We consider a generic particle trapped in an effective harmonic potential. For simplicity we shall consider the 
one-dimensional case, even if the method can be in principle generalized to the three-dimensional case. This particle 
can be an ion trapped by a rf-trap Q or a neutral atom in an optical trap pd] , p^ |. Our scheme however does not 
depend on the specific trapping method employed and therefore we shall always refer from now on to a generic trapped 
"atom" . 

The trapped atom of mass m, oscillating with frequency uj a along the x direction and with position operator 
x = xo(a + at), xq = (?i/2mw a ) 1 / 2 , is coupled to a standing wave in a cavity with frequency ujb, wave-vector k along 
x and annihilation operator b. The standing wave is quasi-resonant with the transition between two internal atomic 
levels |±) separated by Two®. We consider also an external driving of the standing wave with a laser at frequency u>b 
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and of the atomic center-of-mass motion with a classical electric field along the x direction, with frequency lja- The 
resulting Hamiltonian of the system is 

H = ~^~ a z + fiuo a o) a + hujbb^b + ihe(<j+ + &-){b — b^) sin (kx + 4>) 

~qEx Q (a + a f ) sm(u A t + 6) + ih (Be^^tf - B*e lUBt b) , (1) 

where a z = |+)(+| — |— )(— |, <J± = |±)( z f| j an d c is the coupling constant. 

In the interaction representation with respect to Hq = Tllub (b'b + %-J , and making the rotating wave approximation, 
that is, neglecting terms rapidly oscillating at the driving laser frequency uib, this Hamiltonian becomes 

H = a z + huo a a) a + U(u>t, — u>s)b^b + ihe(a + b — er_&t) sin (kx + cf>) 

-qEx Q (a + a [ )sm(u) A t + 6) + ih(Bb i - B*b) , (2) 

where A = ujq — ujb is the atomic detuning. This detuning can be set to be much larger than all the other parameters 
A 3> e, Ldb — u>b, and in this case, the excited level can be adiabatically eliminated, so to get the following effective 
Hamiltonian for the vibrational motion of the atom and the standing wave mode alone |l3|] 

H = H (ujt, — lob) b^b + Tioj a a)a — h— b^b sin 2 (kx + 4>) 

-qEx (a + a r )sm(uj A t + d) + ih(Btf -B*b) . (3) 

If we set the spatial phase = 0, and assume the Lamb-Dicke regime, one can approximate sin 2 (kx + <$>) ~ k 2 x 2 in 
Eq. (|^). Then, in the interaction representation with respect to Tiuj A a) a and making the rotating wave approximation, 
i.e., neglecting all the terms oscillating at u>a (which is of the order of 1 Mhz) or faster because we are interested in 
the dynamics at much larger times, we finally get 

H = H (ujb — ujq — G/2) b^b + fr(uj a — ti>^)a^a — hGrfba^a + ih (Aa^ - A* a) + ih (Btf - B*b) , (4) 

where G — 2(ekxo) 2 / A and A = — qEe~ t6 /2h. This Hamiltonian gives rise to a crossed Kerr-likc effect which could be 
exploited to generate nonclassical states analogously to the all-optical case proposed in Ref. . A similar approach 
was used for Schrodinger cat motional states of atoms in cavity QED |lfj . Here, instead, we are looking for stationary 
nonclassical states. The evolution of the density matrix D of the whole system (vibrational degree of freedom plus 
the cavity mode) is determined by the Hamiltonian (^) and by the terms describing the photon leakage out of the 
cavity with decay rate k and the coupling of the vibrational motion with the thermal environment CthD, that is 

D = L th D- ^[H,D] + ^(2bDtf -tfbD- Dtfb) . (5) 

For the determination of the damping term of the vibrational motion CthD, we note that it occurs at a frequency u> a of 
the order of MHz and that the corresponding damping rate 7 is usually much smaller || . It seems therefore reasonable 
to use the rotating wave approximation in the interaction between the atom center-of-mass and its reservoir, leading 
us to describe the damping of the vibrational degree of freedom in terms of the quantum optical master equation (at 
nonzero temperature) [ p0[ , 

CthP = ^ fa + 1) (2apa^ — a) ap — pa* a) + —n (2a) pa — aa^ p — paa 1 ) , (6) 

where n — [exp (hv/ksT) — 1] 1 , is the number of thermal phonons (ks is the Boltzmann constant and T the 
equilibrium temperature). An analogous treatment is considered in ||. We have to remark, however, that the 
damping and heating mechanisms of a trapped atom are not yet well understood || and that different kinds of 
ion-reservoir interaction have been proposed plf . 

The quantum Langevin equations pof corresponding to the master equation (^|) reads 

b = -i(uj b -ujb- G/2)b + iGa^ab - |fe + B + y^b m (t) , (7) 

T 

a = —i(uj a — uj A )a + iGtfba — —a + A + y^a in (t) , (8) 
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where the input quantum noises bi n (t) and a jn (i) have zero mean and the following correlation functions 

(M*)M*0> = (*»L(*)M*0} = o (9) 

(M*)4(*')> =*(*-*') (10) 
(a m (t)a m (t')> =0 (11) 

(aj n (t)<z* n (t')> (12) 

(a in (t)a] tl (*')> = (» + !)*(*-*')• (13) 

When the external driving terms described by A and B are sufficiently large, the stationary state of the system is 
quasi-classical, that is, the standing wave is approximately in a coherent state with a large amplitude f3 3> 1, and the 
atomic vibrational motion along x is approximately in a coherent state with a large amplitude a>l. The values of a 
and (3 are given by the solutions of the coupled nonlinear equations given by the semiclassical version of the quantum 
Langevin equations Eqs. (Jfy and (||): 

= -i {co b -ujB- G/2 - G\a\ 2 ) + (14) 

= -i(iv a -u; A -G\P\ 2 )a~ ^a + A. (15) 

Since a ~ 2A/"f and j3 ~ 2B/k, the semiclassical condition for the steady state is satisfied when .4^7 and B k. 

The fluctuations around this steady state are instead described by quantum mechanics and their dynamics can be 
obtained by appropriately shifting both modes, i.e., b — > b + (3 and a — > a + a. In the semiclassical limit \a\, \0\ ^> 1 
it is reasonable to linearize the equations, and since it is always possible to tune the two driving frequencies loa and 
ujb so to have zero detunings, i.e., \a\ 2 — (ujb — lob)/G — 1/2, \(3\ 2 — (u> a — u>a)/G, the linearized quantum Langevin 
equations for the quantum fluctuations around the steady state can be written as 

b = iGf3{a*a + aa))-^b + ^b m {t), (16) 

a = iGa{p*b + 0)- ~a + ^a in (t). (17) 

The effective linearized Hamiltonian leading to Eqs.(|l6|), (|l7|), can be written as 

H = h x YX , (18) 

where \ = — 4G|a| |/5| , Y = (6e~ 1 ^ +b^e i ^)/2 is the standing wave field quadrature with phase <f)p equal to the phase 
of the classical amplitude (3, and X — (ae _i< ^ Q + a' ! e 1 ^ a )/2 is the vibrational quadrature with phase <j> a given by the 
phase of the classical amplitude a. For the sake of simplicity we shall consider <p a — 4>p ~ 0, i.e., the atomic position 
quadrature, from now on, even if the following considerations can be easily extended to the case of generic phases. 
Note that in order to remain in the Lamb-Dicke regime it is required that fcxo|a| <^ 1; however the linearisation is 
justified only when \a\ ^> 1, and therefore we need kxo <C -Jj <S 1. 

Eq. ( |l8|) implies that an effective continuous, quantum non-demolition (QND) measurement of the phonon quadra- 
ture X is provided by the homodyne measurement of the light outgoing from the cavity, which plays the role of the 
"meter". In fact, the homodyne photocurrent is jicf 

I(t) = 2 V K(Y v (t)) e + y/TjKS(t), (19) 

where Y v = (be~ lip + b^e lv )/2 is the measured quadrature, the phase ip is related to the local oscillator, and r\ is the 
detection efficiency. The subscript c in Eq. ( |l9| ) denotes the fact that the average is performed on the state conditioned 
on the results of the previous measurements and £(£) is a Gaussian white noise |l6| ]. In fact, the continuous monitoring 
of the field mode performed through the homodyne measurement, modifies the time evolution of the whole system. 
The state conditioned on the result of measurement, described by a stochastic conditioned density matrix D c , evolves 
according to the following stochastic differential equation (considered in the Ito sense) 

D c = C th D c - ~ [H, D e ] + I (2bD c tf - tfbD c - D c tfb) 

+ y/VK&t) (e-^bDa + e^Dctf - 2{Y V ) C D C ) . (20) 
We note that by performing the average over the white noise £ (i) , one gets the master equation of Eq. (H) . 
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It is now reasonable to assume that the standing wave mode is highly damped, i.e. k ^> \ (this does not conflict 
with the preceding assumptions, since the coupling constant \ = —8e 2 (kx ) 2 \af3\/ A is usually smaller than the cavity 
decay rate). This means that the radiation field will almost always be in its lower state |0)fc (displaced by an amount 
(3). This allows us to adiabatically eliminate the field and to perform a perturbative calculation in the small parameter 
obtaining (see also Ref. |l7]]) the following expansion for the total conditioned density matrix D c 

D c = (p c - ^X Pc x) <g> |0) 66 (0| - (Xp c <g> |1) M (0| - Pc X ® |0) 66 (1|) 

\ K z J K 

2 2 

+ \Xp c X ® |1) 66 (1| - (X 2 Pc (8 |2) 66 (0| + p c X 2 ® |0) w <2|) , (21) 



where p = Trj D is the reduced density matrix for the vibrational motion. In the adiabatic regime, the internal 
dynamics instantaneously follows the vibrational one and therefore one gets information on X by observing the 
quantity Y v . The relationship between the conditioned mean values follows from Eq. ( |2l| ) 

{Y v {t)) c = ^{X(t)) c sm<p. (22) 

Moreover, if we adopt the perturbative expansion ( pl| ) for D c in ( pfj| ) and perform the trace over the internal mode, 
we get an equation for the reduced density matrix p c conditioned to the result of the measurement of the observable 
(Y v (t)) c , and therefore (X(t)) c 

x 2 



Pc = C thPc - ±- [X, [X,p c ]} + y/vx 2 /*^) (ie iv p c X - ie-^Xpo + 2sin^(X(t)) c p c ) . (23) 
^k 

This equation describes the stochastic evolution of the vibrational state of the trapped atom conditioned to the result 
of the continuous homodyne measurement of the light field. The double commutator with X is typical of QND 
measurements. 

III. THE FEEDBACK LOOP 

We are now able to use the continuous record of the atom phonon quadrature to control its motion through the 
application of a feedback loop. We shall use the continous feedback theory proposed by Wiseman and Milburn |lS| ]. 

One has to take part of the stochastic output homodyne photocurrent I(t), obtained from the continuous monitoring 
of the meter mode, and feed it back to the vibrational dynamics (for example as a driving term) in order to modify 
the evolution of the mode a. To be more specific, the presence of feedback modifies the evolution of the conditioned 
state p c (t)- It is reasonable to assume that the feedback effect can be described by an additional term in the master 
equation, linear in the photocurrent I(t), i.e. 

\Pc(t)] fb = ^ Z ^JCpS) t (24) 

3 vx 

where r is the time delay in the feedback loop and K. is a Liouville superoperator describing the way in which the 
feedback signal acts on the system of interest. 

The feedback term (Q) has to be considered in the Stratonovich sense, since Eq. (^4]) is introduced as limit of a 
real process, then it should be transformed in the Ito sense and added to the evolution equation (|23|). A successive 
average over the white noise £(t) yields the master equation for the reduced density matrix p = TibD in the presence 
of feedback. In the general case of a nonzero feedback delay time, one gets a non-Markovian master equation which is 
very difficult to solve jl8j (see however Ref. [^9|). Most often however, the feedback delay time is much shorter than 
the characteristic time of the a mode, which in the present case is given by the energy relaxation time 7 -1 , and in 
this case the dynamics in the presence of feedback can be described by a Markovian master equation |18| ], which is 
given by 

p = C thP - J- [A, [A, p]\ + K (ie**pX - ie-*Xp) + t^ttP- (25) 

The third term is the feedback term itself and the fourth term is a diffusion-like term, which is an unavoidable 
consequence of the noise introduced by the feedback itself. 
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Then, since the Liouville superoperator JC can only be of Hamiltonian form we choose it as JCp = g [a — a^ , p] /2 
|l7| , which means feeding back the measured homodyne photocurrent to the vibrational oscillator with a driving term 
in the Hamiltonian involving the quadrature orthogonal to the measured one; g is the feedback gain related to 
the practical way of realizing the loop. One could have chosen to feed the system with a generic phase-dependent 
quadrature, due to the homodyne current, however, it will turn out that the above choice gives the best and simplest 
result. Since the measured quadrature of the vibrational mode is its position, the feedback will act as a driving for 



the momentum. Using the above expressions in Eq. (25) and rearranging the terms in an appropriate way, we finally 
get the following master equation: 



where 
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Eq. (|26j) is very instructive because it clearly shows the effects of the feedback loop on the vibrational mode a. 
The proposed feedback mechanism, indeed, not only introduces a parametric driving term proportional to gshiip, 
but it also simulates the presence of a squeezed bath, characterized by an effective damping constant T and by the 
coefficients M and N, which are given in terms of the feedback parameters [0. An interesting aspect of the effective 
bath described by the first four terms in the right hand side of (|2^) is that it is characterized by phase-sensitive 
fluctuations, depending upon the experimentally adjustable phase 09. 



IV. THE STATIONARY SOLUTION 



Because of its linearity, the solution of Eq. (|26j) can be easily obtained by using the normally ordered characteristic 
function [Eof C(X,\*,t). The partial differential equation corresponding to Eq. (p6[) is 



(30) 



dt + ^Xd x + ^X*d x , +|sin^(A9 A . + A*d A ) j C(A, A*, i) 

-rW|A| 2 + ^M+|sin^ (A*) 2 + ^Ar + |sin^ X 2 } C< A. A : . / ) . 



The stationary state is reached only if the parameters satisfy the stability condition, i.e. <?sino9 < 7. In this case 
the stationary solution has the following form 



C(A, A*, 00) = exp 



-C|A| 2 + ^(A*) 2 + ^*A 2 



where 



c = 



AT 2 + 3 sino9(rRe{M} + 2vlm{M}) + g 2 sin 2 ip/2 
T 2 — g 2 sin 2 09 



(N + l/2) 5 sino9 + TRe{Af} (r 2 - g 2 sin 2 09)Im{M} 
p = r — — - ; n h 1- 



r 2 — a 2 sin 2 



r 2 — a 2 sin 2 



'•P 



(31) 

(32) 
(33) 



Under the stability conditions and in the long time limit (t — > 00) the variance of the generic quadrature operator 
X e = {ae lt> + a)e- l0 )/2 becomes 



(XI) = 



( + Re{pe 2tti } 



(34) 
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V. DISCUSSIONS AND CONCLUSIONS 



For the X quadrature Eq.(|3J) can be simply written as 

(35) 



1 

2 + n eJ7 



with n e ff = C + Rc{/i}. In absence of feedback (g = 0) we have n e ff = n, otherwise n e ff can be smaller than n, 
providing a stochastic localisation in the X quadrature. Depending on the external parameters, it can also be negative 
(but it is always n e ff > — 1/2) accounting for the possibility of going beyond the standard quantum limit. This is a 
relevant result of the present feedback scheme since it is able to reduce not only the thermal fluctuations but even 
the quantum ones. 

The potentiality of this feedback mechanism is clearly shown in Fig.l, where n e ff goes well below zero for increasing 
values of x- Instead in Fig. 2 we have sketched the phase space uncertainty contours obtained by cutting the Wigner 
function corresponding to Eq.(|T]) at l/y/e times its maximum height. We see that the state resulting from the 
feedback action (solid line) has a relevant contraction in the Xo=o direction, but the same uncertainty in the Xg =1v /2 
direction with respect to the state of the system undergoing measurement without feedback (dashed line). We refer 
to this type of noise reduction produced by feedback loop as squashing, whereas squeezing refers to conventional 
quantum noise reduction ]l0| . 

Summarizing, we have proposed a feedback scheme based on an indirect (QND) measurement which is able not only 
to contain the heating of the vibrational motion of a trapped ion, but also to produce nonclassical motional states 
(squashed ones) . Up to this point we have not discussed the specific way in which a particular feedback Hamiltonian 
could be implemented. In our case, it is important to be able to realize a term in the feedback Hamiltonian proportional 
to the quadrature orthogonal to X. This is not straightforward, but could be realized by using the feedback current 
to vary an external potential applied to the atom without altering the trapping potential p2| ] . In principle the model 
could be extended to the three dimensional case, one should only consider three orthogonal standing waves far from 
resonant transitions. 

In conclusion, although the experimental implementation of the presented model may not be easy, it is certainly 
a promising experimental challenge, stimulated by the possibility of producing nonclassical states for trapped atoms 
and of controlling their heating to minimize decoherence effects, especially in quantum information processing f23|. 
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FIG. 1. The quantity n e // is plotted vs g for different values of x\ the values of the other parameters are n = 0.5, 7 = 10 2 
s _1 , k — 10 2 s" 1 and r\ = 0.8. The quantities \ an d 9 are expressed in s _1 . 
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Xe = 7l/2 
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FIG. 2. The phase space uncertainty contours are represented for g = (dashed line) and g — 0.025 s 1 (solid line), the 
values of the other parameters are n = 0.5, x = 2.5 s-\ 7 

= 10 2 s 1 , k = 10 2 s 1 and r\ = 0.8. The dotted line represents 

the vacuum noise. 
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